Spin-Peierls transition in NaV205 in high magnetic fields 
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We investigate the magnetic field dependence of the spin-Peierls transition in NaV2 05 in the field 
range 16T-30T. The transition temperature exhibits a very weak variation with the field, suggesting 
a novel mechanism for the formation of the spin-Peierls state. We argue that a charge ordering 
transition accompanied by singlet formation is consistent with our observations. 



PACS: 75.30.Kz, 75.40.Cx, 75.50.Ee 

The Peierls instability takes place in one-dimensional 
systems and can give rise to complex and fascinating be- 
havior. In itinerant electronic systems the instability is 
driven by the coupling of electrons to the phonons of the 
lattice . Any coupling at T = leads to the formation 
of the Peierls state which is characterized by charge or- 
dering (gap in the electronic spectrum) and a finite lattice 
distortion. Similar phenomena occur in purely insulating 
spin systems, where the spin-phonon coupling is responsi- 
ble for the formation of a singlet ground state with neigh- 
boring spins pairwise bound into singlets Q. The spin- 
Peierls ground state shows a characteristic gap in the ex- 
citation spectrum and has been observed in a variety of 
organic compounds, such as (TTF)[CuS4C4(CF3)4] §• 
At high temperatures these materials behave as non- 
interacting Heisenberg chains, while below the transi- 
tion temperature, T^, the magnetic exchange acquires 
an alternating component. In 1993 the first inorganic 
spin-Peierls compound CuGeOa was discovered with 
Tc ~ lAK . This material, like its organic predecessors, 
shows a characteristic Id Heisenberg (Bonner-Fisher)-like 
magnetization at high T with a sharp drop at Tc, indicat- 
ing a non-magnetic ground state. Very recently, a second 
inorganic compound, NaV205 was shown to behave as a 
spin-Peierls material with Tc « 34K Q]. The properties 
of NaV205 however have proven to be quite controver- 
sial, thus stimulating the research reported in this letter. 

Magnetic susceptibility measurements of NaV205 in- 
dicate a transition to a non-magnetic phase at Tc [§,1). 
This can be understood within the framework of a spin- 
phonon coupling driven transition on a Heisenberg chain 
The antiferromagnetic exchange, J, was estimated 
to be Jw560K. The low-temperature structure which is 
assumed in this interpretation of the data, is that of mag- 
netic chains formed by the spin 1/2 ions along the 
crystalline b-axis, separated by spinless chains. This 
scenario implies a lattice distortion in one direction only. 
However, recent experiments have shown that the above 
picture is not satisfactory. X-ray diffraction measure- 



ments indicated that the system should be viewed as a 
quarter-filled ladder made of V*'^^ chains ^,@|, meaning 
that a spin of 1/2 is not attached to a single V ion, but 
rather to a rung of the ladder, i.e. a V-O-V orbital. Sub- 
sequent NMR j|] analysis revealed that below Tc, two 
inequivalent types of V sites - V''+ and V^"*" appear, sug- 
gesting that charge ordering occurs in the spin-Peierls 
phase. Charge disproportionation leaves room for period 
doubling in more that one crystallographic direction, con- 
sistent with additional X-ray (9|jl^ and NMR studies. 
These works suggest that lattice distortion takes place 
in the (a,b) plane (where b is the direction along the 
chains and a is perpendicular to the chains). A number 
of theoretical studies p^-p7[| have addressed the possibil- 
ity of charge ordering in 1/4 filled systems, where both 
electron-lattice and electron-electron interactions are in- 
cluded. The most probable scenario at present seems to 
be the "zig-zag" order proposed in Ref. |l2j where the 
charge density (i.e. the sites V'''^** with deviation 5 
from the average valence) is distributed in a zig-zag fash- 
ion along the ladder direction. As emphasized in Ref. p^ ] 
the Coulomb repulsion in combination with the electron- 
lattice interaction can drive such a transition, while the 
formation of a spin singlet ground state "follows" the 
charge order. Charge modulation is consistent with the 
analysis of the observed magnetic excitation spectra , 
Raman spectra p9(|, a s well as the anomalies in the ther- 
mal conductivity |2(| and the dielectric constant at 
Tc. 

The present work attempts to gain further insight into 
the nature of the spin-Peierls transition in NaV2 05 by 
addressing the magnetic field dependence of the transi- 
tion temperature in very high fields. Previous studies in 
fields up to 5.5T ||] have found behavior consistent with 
the theoretical predictions and similar to the previously 
known spin-Peierls compounds [Q. However, subsequent 
measurements in higher fields, up to 14T [ p^ and 16T 
p3[ , have found much weaker field dependence. These 
experiments were based on a determination of Tc from 
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the changes in the elastic constants [g2| and the specific 
heat unHke the measurement in Ref. Q which de- 
termined Tc from the drop in the magnetization. 

In this work we have measured the magnetization of 
two NaV205 single crystals in magnetic fields from 16T 
to 30 T. The crystals were grown by high temperature so- 
lution growth from a vanadate mixture flux. The masses 
of the samples investigated were 1.9 and 3.1 mg respec- 
tively and they had irregular parallelepiped shapes with 
smooth, faceted faces. The single crystals were charac- 
terized with an Enraf-Nonius CAD4 single crystal diffrac- 
tometer using Mo-radiation. The results of the structure 
refinement were the same as reported earlier in Ref. [Q. 
Magnetization was measured using a standard metal foil 
cantilever beam magnetometer. The "T" shaped flexi- 
ble cantilever beam was made from a 7.62 /im thick heat 
treated MP35N alloy. The dimensions of the "T" were 
approximately 8 mm on a side. The gap between the "T" 
and the parallel fixed reference electrode was approxi- 
mately 800 /im. The sample was mounted using a small 
amount of vacuum grease. In the presence of a DC mag- 
netic field the interaction of the magnetic moment of the 
sample with the field results in a force and/or torque, de- 
flecting the beam and changing the capacitance between 
electrodes. A capacitance bridge was used to monitor the 
changes in force (magnetization) for temperature sweeps 
in fixed field. Since MP35N is magnetic (typically 13.5 
/iemu/g at 78 K), the same bare cantilever was measured 
under the identical conditions (sweep direction and sweep 
rate) as the cantilever -|-sample combination to provide a 
background reference. The temperature dependence of 
the cantilever capacitance was compensated for in the 
same way. A room temperature measurement of the can- 
tilever's sensitivity showed that a force of 3 nN could be 
resolved. 

Cantilever displacement can arise from either a torque 
or a force on a sample with a magnetic moment. When 
the sample is at field center, where the field gradient 
is zero, then torque (c>c m x B) will dominate. Strictly 
speaking, if the sample is isotropic and there are no shape 
factors, then there is no torque on the cantilever for fields 
applied along the direction of displacement (perpendicu- 
lar to the sample). On the other hand, when the sample 
is raised (or lowered) away from field center, the force 
term (F — m dB/dz oc x B dB/dz) will usually dom- 
inate, although torques can still be present. Figure ^ 
shows temperature sweeps taken for the 1.9 mg sample 
at the three indicated fields in the legend. The cantilever 
was located in a position where the field gradient was 
maximum. The maximum field at this position (24T) is 
80% of the field center maximum (30T). The change in 
capacitance, AC, which is proportional to the change in 
magnetization, is calculated at each field by subtracting 
the background trace (cantilever alone) from the sample 
trace (sample-f-cantilever) . This quantity is divided by 
and plotted as the ordinate in Fig. As seen from 



the figure, the data scale reasonably well for the three 
different fields, confirming the dependence expected 
from both torque and force contributions. To accentuate 
the small shift in the transition temperature, we plot in 
the inset the derivative of AC/S^ with respect to the 
temperature. From the position of the peaks we can de- 
termine the field-dependent transition temperature. 

To measure Tc at the maximum field of 30T, the sam- 
ple was placed at field center (sensitive to the torque 
only) and the data collected and analyzed as described 
above. A similar scaling with B^ was observed. Figure ^ 
shows the derivative of AC with respect to temperature 
at field center. Plotted in this way, the shift in Tc can 
be clearly seen. Data similar to those plotted in Figs. 
1 and 2 were obtained for a second sample with mass 
3.1 mg, and for reversed fields. In all cases the shifts 
in Tc were equal to or less than the shifts shown in Fig. 
^. In Fig. H (inset) we also show the results of magne- 
tization measurements at low fields using a commercial 
SQUID magnetometer (MPMS7). The singlet formation 
at Tc is clearly observable, but no shift of Tc can be ob- 
served within measurement accuracy in fields up to 5T, 
in agreement with previous work 

In Fig. ^ we present our high field data for the vari- 
ation of Tc in terms of ATc/Tc(0) = Tc(i/)/Tc(0) - 
1 and the square of the scaled magnetic field h — 
gpLBH/2kTc{Q) iQ. This scafing is expected in spin- 
Peierls systems, and for small fields <C 1, the relative 
variation of Tc should be quadratic Q : 

ATc/Tc(0) = -ah^. (1) 

The data of Fig. ^ follows this dependence quite well, 
and we estimate aexp ~ 0.072(8). The value Tc(0) was 
not measured directly but was estimated from an extrap- 
olation to zero field of the quadratic dependence of Tc 
vs. H to be Tc(0) = M.2K. This value is close to 
published values and to the Tc(0) measured by us us- 
ing SQUID magnetometer measurements of the magne- 
tization of a 40mg polycrystalline sample. The combi- 
nation of our high field data and the lower field data 
of previous measurements gives the variation of Tc over 
a large range of magnetic field and shows a very weak 
dependence. In contrast, the "conventional" inorganic 
spin-Peierls compound CuGeOs exhibits a much stronger 
field dependence with a = 0.39 in good agreement 
with the theory. The theoretical values of as p predicted 
for the spin-Peierls transition are asp = 0.44 or 0.36, 
depending on the way interaction effects are taken into 
account The first, larger number corresponds to the 
Hartree approximation for the interactions between the 
Jordan- Wigner fermions, representing the localized spins 
p6[ . The value 0.36 is obtained by exact treatment of 
the correlation effects which is possible in the Lut- 
tinger liquid framework in one dimension p8[ . In both 
cases the characteristic scaling H/Tc{0) which appears in 
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Eq.(l) is due to the commensurate nature of the dimer- 
ized phase. For large fields, corresponding to a reduction 
of Tc by a factor of TclTc{Q) — 0.77, a transition into 
an incommensurate phase is expected to take place [ p7[ . 
Such a transition is less sensitive to magnetic field and 
has been observed in a variety of spin-Peierls materials 
1^,^. In NaV2 05 however, a transition into such a mod- 
ulated phase does not seem to take place, since even in 
the highest field (SOT), rc(30T)/Tc(0) = 0.97, which is 
very far from the expected incommensurate boundary. 
Notice that even in a field as high as SOT the scaled ratio 
h — 0.59 is quite small due to the large Tc(0). 

We now discuss the possible sources for the difference 
between the measured value aexp and the theoretically 
predicted one asp = 0.S6 ~ 5aexp for spin-Peierls sys- 
tems. In addition to this discrepancy, any theory of 
NaV2 05 should also be able to explain the large value 
of the ratio 2A/rc(0) « 6 (A « lOOK @ being the spin 
gap), where a mean-field value of S.52 might be expected. 

As discussed in the introduction, a transition into a 
charge ordered state in a 1/4- filled system is consistent 
with a number of recent experiments. Although it is not 
clear whether the charge density wave (CDW) precedes 
or forms simultaneously with the magnetically dimer- 
ized spin-Peierls state, it seems certain that the physics 
of charge ordering must be taken into account. Recent 
numerical work has shown |l6|Jl7| ] that CDW and spin- 
Peierls order can co-exist in quasi one-dimensional 1/4- 
fiUed electronic systems. If we assume that the CDW 
formation is the driving force behind the opening of a 
spin gap, as argued in Ref. |T^, then the "charge" part 
of the transition will be mainly responsible for the Tc{H) 
dependence. In a system of non-interacting electrons, 
undergoing a Peierls transition into a (commensurate) 
CDW state, the decrease of Tc for small magnetic field 
(coupled to the electron spin via a Zeeman term) is also 
described by Eq.(l), but with acDW = 0.21 

Two effects, orbital coupling and electron-electron in- 
teractions, could further modify this result. Orbital ef- 
fects are known to be present when nesting is imperfect, 
and generally compete with the Pauli terms, producing 
a flatter dependence of Tc on H, i.e. a further reduction 
of acDW |l30|| . However spin-orbit interactions lead to 
anisotropic variation of Tc with respect to the magnetic 
field direction. In NaV2 05 this variation has been found 
to be extremely weak P,p|, p2|j23| , which is also confirmed 
in this work, and consequently the orbital effects can be 
ruled out as a source of the weak Tc{H) dependence. On 
the other hand, electron-electron interaction effects do 
not reflect anisotropics and are important in the forma- 
tion and stabilization of a CDW state Jl^ , p^ . In gen- 
eral, the stability of the CDW depends on the strength 
of the electron-phonon coupling (which drives the tran- 
sition) and on the on-site and nearest-neighbor Coulomb 
correlations [ p^ . 

To demonstrate this latter point concerning strong cor- 



relation effects, we consider the simplified model of a 
Hubbard chain with an on-site repulsion U. We treat 
the phonons adiabatically, as in Ref. but take into 
account the electron-electron interaction following Ref. 
p8[ , i.e. calculate the polarization bubble exactly for the 
Luttinger liquid. In this case it is known that Tc{H ~ 0) 
increases with respect to its value aiU = Q . For finite 
magnetic field we find, at J7 ~ 2t (where t is the band- 
width), that the coefficient a drops to acDW ~ 0.15, 
i.e. below the non-interaction value of 0.21. This is not 
surprising and in fact is quite similar to the difference 
between the mean-field and the exact treatment in the 
spin-Peierls case {asp = 0.44, 0.36, respectively). The 
essence of the effect is in the different type of diver- 
gence in the polarization bubble with and without inter- 
actions. While in the free case the polarization diverges 
logarithmically at small frequencies, in a Luttinger liq- 
uid the stronger, power law dependence sets in ||2^, and 
the Peierls instability is effectively enhanced. Thus the 
interaction effects, being naturally more important for 
the CDW formation (compared to the spin-Peierls case), 
can produce a weaker Tc{H) dependence. A more real- 
istic calculation based on a Hamiltonian appropriate for 
NaV205 would be very desirable. 

The orbital and interaction effects discussed above are, 
strictly speaking, valid only for an isolated chain. It 
was assumed that inter-chain interactions are sufficiently 
strong to suppress the fluctuation effects, typically im- 
portant in one-dimensional systems p^ . The fluctua- 
tions arc known to reduce T'c(O) below the mean- field 
value and cause a specific heat jump at the transition 
Acp several times the mean-field one. The large ob- 
served ratio A/rc(0) (twice the mean- field), in combi- 
nation with a Acp about ten times the mean-field value 
1^ suggest that fluctuations indeed could be important 
in NaV2 05. At the same time one should have in mind 
that, due to the specific structure of NaV205, transverse 
interchain interactions are expected to play a crucial role 
in the stabilization of the ordered phase, in particular 
the formation of the spin gap and doubling of the period 
in the (a,b) plane |9|jri[]. The vanadium displacements 
are nearly absent along the ladder direction (b-axis), and 
largest perpendicular to the ladder direction both along 
the a- and c-axis [|lO|. Thus it is not clear whether fluc- 
tuation effects have to be necessarily invoked to explain 
the large A/Tc(0) ratio in this material as is traditionally 
done, or whether the large A/Tc(0) ratio is intimately re- 
lated to the anomalously weak variation of Tc with field 
reported in this work. 
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FIG. 1. Change in capacitance measured with a cantilever 
beam magnetometer off field center for a NaV205 single crys- 
tal in B = 16T, 20T, and 24T. AC is proportional to the 
magnetization of the sample and has been normalized by the 
square of the magnetic field. The inset shows the derivative of 
the scaled, background subtracted data with respect to tem- 
perature. Tc is determined from the position of the peaks. 



FIG. 2. Derivative of the unsealed capacitance readings 
(proportional to magnetization) with respect to temperature 
for the sample shown in Fig. 1 located at field center. Inset 
shows low field SQUID magnetization measurements. 



FIG. 3. Relative variation of Tc as a function of the scaled 
magnetic field h = gfisH /2kTc{0) (see text). The circles 
are our data (numbers represent the values of the field in 
Tesla), and the squares are data from Ref. based on 

measurements of the specific heat jump at the transition. 



4 



1.4x10 - 

1.2x10"^- 

1.0x10 ^ - 

8.0x10'- 

6.0x10"'- 

4.0x10"'- 

2.0x10"'- 
0.0 - 

-2.0x10"'- 




16T 
20T 
24T 





5.0x10"' - 




4.0x10"' - 


1- 






3.0x10"' - 


CM 

m 




o 


2.0x10"' - 


<, 










I.Oxlo' - 




0.0 - 



1.0x10 




24 26 28 30 32 34 36 38 40 42 
T(K) 



— ' — I — ' — I — ' — r 

24 26 28 30 



T 



32 34 

T(K) 



1 — ' — I — ' — I — ' — I — ' — 

36 38 40 42 



8.0x10" 



6.0x10 - 



^ 4.0x1 0"'- 

LL 
Q_ 



Q 2.0x1 0"H 

<l 



0.0- 



-2.0x10" 



W^ JV" •••• /♦ 




20 T 
25 T 
SOT 



1T 
5T 



» / 



24 26 28 30 32 34 36 38 40 42 
T(K) 



— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T" 

24 26 28 30 32 34 36 38 40 

T(K) 



42 



